Elastic properties of 2D colloidal crystals from video microscopy 
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Elastic constants of two-dimensional (2D) colloidal crystals are determined by measuring strain fluc- 
tuations induced by Brownian motion of particles. Paramagnetic colloids confined to an air-water 
interface of a pending drop are crystallized under the action of a magnetic field, which is applied 
perpendicular to the 2D layer. Using video-microscopy and digital image-processing we measure 
fluctuations of the microscopic strain obtained from random displacements of the colloidal particles 
from their mean (reference) positions. From these we calculate system-size dependent elastic con- 
stants, which are extrapolated using finite-size scaling to obtain their values in the thermodynamic 
limit. The data are found to agree rather well with zero-temperature calculations. 
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During the last two decades interest in colloidal sys- 
tems has grown substantially, on one hand because of 
their widespread technological applications and on the 
other due to the availability of precisely calibrated par- 
ticles for use as model systems for studying phenomena 
in classical condensed matter physics [1]. The crystal- 
lization of colloids, both in two and three dimensions 
has been a continuous matter of interest. The research 
mostly focused on the analysis of structure and dynamics 
of colloidal systems on different length and time scales 
through static or dynamic light scattering techniques. 
Measurements of elastic constants of colloidal crystals, 
however, have been limited to the determination of the 
shear modulus [i. This was based on the observation of 
shear induced resonance of a crystal using light scattering 
techniques (see [2] for a recent work) . The value of \x is 
found to depend strongly upon the crystalline morphol- 
ogy and changes significantly between randomly oriented 
crystallites and shear-ordered samples [3]. In addition, 
using this method only a very reduced number of modes 
can be investigated. Very recently, the elastic moduli of 
colloidal solids have also been estimated [4] by observ- 
ing relaxation behaviour after deformations using laser 
tweezers. 
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FIG. 1. A snapshot of the triangular lattice of paramag- 
netic colloidal particles. A few thousand snapshots such as 
this, taken at regular time intervals of about one second were 
used to calculate elastic constants. 

In this Letter, we report an experimental determination 



of the equilibrium elastic properties of two -dimensional 
(2D) colloidal crystals from "snapshots" of particle po- 
sitions obtained using video microscopy. The present 
method is completely non-invasive, accurate and free 
from any adjustable parameters. 
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FIG. 2. The comparison of the measured elastic con- 
stants (in units of kT/a 2 , a is the lattice parameter) to the 
zero-temperature calculation (solid lines) reveals an very good 
agreement. Note that the shear modulus fj, is multiplied by 5 
for reasons of clarity. 

The mechanical properties of a macroscopic solid, ac- 
cording to classical elasticity, are well described by a 
small set of elastic constants. These can be measured 
by the strain-response of the solid under the application 
of an appropriate (macroscopic) stress. On a mesoscopic 
scale, which is still sufficiently coarse-grained to apply 
elasticity theory but small enough such that the Brow- 
nian motion of the particles is observable, thermally in- 
duced strain fluctuations can be used to determine elastic 
constants. We have carried out a detailed study of these 
strain fluctuations in a two -dimensional (2D) colloidal 
crystal by recording the microscopic positions of the par- 
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tides within a square cell (of size L) containing a defect - 
free single crystal (see Fig. 1). Recordings were made 
at regular time intervals — large compared to typical 
correlation times of the colloid. Uncorrelated snapshots 
obtained are analyzed to calculate the average particle 
positions — the "reference" lattice. The gradients (ob- 
tained by finite differences) of the displacement vectors 
then yield the microscopic strains. These microscopic 
strains are used to obtain strain fluctuations over a hier- 
archy of length scales corresponding to smaller sw&-cells 
of size L b < L contained within our cell. The width of 
the probability distributions of the strains are related to 
the elastic constants Cijki (Lb) obtained as a function of 
Lb which may subsequently be extrapolated, using a sys- 
tematic finite size scaling analysis [5]. The macroscopic 
(L — > oo ) values of these quantities, thus obtained, are 
compared to theoretical predictions without any fitting 
parameters. The central result, the bulk (K) and shear 
(n) elastic moduli are shown as a function of the inter- 
action strength T for our colloidal system (see below) in 
Fig. 2. 

Our experimental setup [6,7] is composed of su- 
per - paramagnetic spherical colloids [8] of diameter 
d = 4.5 fim and mass density 1.7 kg/dm 3 . They are con- 
fined by gravity to a water/air interface, which is formed 
by a cylindrical drop suspended by surface tension in 
a top-sealed ring. The flatness of the water-air inter- 
face (0 = 8 mm) is controlled within ±1 /zm [6]. For 
weak magnetic fields B applied perpendicular to the in- 
terface the induced magnetic moment M depends lin- 
early on B, i.e. M = xB with an effective magnetic sus- 
ceptibility x [6]- The repulsive magnetic dipole-dipole 
potential, between particles i and j separated by a dis- 
tance Tij, V(rij) — Tr^ 3 dominates the interaction 
and is absolutely calibrated by the interaction strength 
T= (/i /47r)(xB) 2 (7rn) 3 / 2 /fcr where n denotes the 2D 
volume-fraction of the particles, k is the Boltzmann con- 
stant, T is the ambient temperature and distances, r^ , 
are in units of the mean interparticle spacing. 

The experiments were carried out as follows: At high 
T, in the crystalline phase, the system was equilibrated by 
the application of small AC magnetic fields in the plane 
of the particles. Eventually a defect-free crystal contain- 
ing several thousand particles is obtained. The entire 
sample consists of approximately 10 5 particles. The co- 
ordinates of typically 1000 particles in the (square) field 
of view were recorded in time. About 500 to 1000 in- 
dependent coordinate sets, each about a second apart, 
are necessary to obtain sufficient statistics. The time in- 
terval corresponds to diffusion over about 1 micrometer 
which is about one pixel — the limit of resolution of our 
digital camera. The results do not vary significantly over 
the range of the number of data sets used. 

After the determination of the mean position of each 
particle R° (taken over the entire set) [9] the instanta- 
neous displacement u(R°) = R R° from the mean po- 
sition was calculated for each frame and for all particles. 



The set of coordinates R° constitutes, therefore, our ref- 
erence lattice and the fluctuating displacement variable 
u(R°) is known at every reference lattice point. Once 
the displacements are known, the elements of the clastic 
strain tensor and the local rotation 9 may be defined 
as gradients of the displacement u over the spatial coor- 
dinates r. 
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When the derivatives are replaced by appropriate finite 
differences, these quantities may be evaluated at every 
reference lattice point R°. In order to evaluate elastic 
constants, the microscopic strains and rotations need to 
be coarse-grained or averaged over a sequence of sub- 
systems (obtained in our case by dividing our square 
cell into integral numbers of smaller square sub-cells) of 
size Lb = L/b, (b = 4,5, .., 18) to obtain size dependent 
strains and rotations e^(Lf,) (and 9b{Lb)) 

e\ i (L b )ML b ) = ± 1 I £ «(r),0(r)dr. (2) 

The fluctuations of these quantities for a particular set of 
sub-systems (all of size Lb) over the different snapshots 
yield elastic constants for the size Lb from standard ther- 
modynamic relations. Although this procedure has been 
used [5] for obtaining elastic moduli of the hard and soft 
disk systems from computer generated configurations, we 
find that some essential modifications are required before 
the methods of Ref. [5] can be applied to our system. 
In a typical computer simulation [5,10-14] the system is 
placed in a hard constraint determined by the ensemble 
used. For example in a constant (zero) strain ensemble 
[5] the displacement field u(r) strictly vanishes at the 
boundary so that the strain, integrated over the entire 
system is strictly zero. Our experimental cell, on the 
other hand, is a small region embedded within a much 
larger crystal. The only constraints which are appropri- 
ate for this case is that the displacement u are continuous 
throughout the system and u(r) — > as r — > oo. In the 
limit of linear elasticity for our system where the elastic 
correlation length is assumed to be much smaller [5] than 
Lb, our problem reduces to the determination of the total 
energy of a single, independent, strain fluctuation in an 
elastic medium subject to these constraints. 

Consider, therefore, a region of size Lb, with a fixed 
constant strain (or rotation) e*-(L&) (6b(Lb)) embedded 
in an (infinite) elastic continuum. This general problem 
has been studied in detail in several standard texts on 
classical theory of elasticity [14,15]. Recall that a 2D 
hexagonal lattice has isotropic elastic behavior [14] and, 
therefore, can be completely described by two indepen- 
dent elastic constants, which we choose to be the bulk 
modulus K and the shear modulus fi. While the for- 
mer is related to the fluctuations of the volume of the 
sub-systems the latter may be evaluated from the fluctu- 
ations of the angle of rotation of the sub-systems. Let us 
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focus on a small sub-cell within our cell. Consider for the 
moment that we have a disk of radius R b for simplicity, 
the final results will be cast in a form independent of the 
shape of the sub-system. Within this disk, the strains 
are given by their values which are the averages over the 
area of the disk. This disk is embedded in an infinite 
elastic medium with elastic moduli K and 

We consider first a homogeneous expansion (or com- 
pression) of the disk by Rb — > Rb + Ar. The correspond- 
ing radial-displacement u r is given as, 



u r = Ar ■ r/Rb, r < Rb 
= Ar-R b /r, r > R b 



(3) 



The angular part u v — by symmetry. The displace- 
ments (u r ,u v ) are related to the strain tensor by the 
following equations [15]: 
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Making use of the (quadratic) free energy density / of 
the elastic continuum, 



f = \\K{e r 



+ e w ) 2 +n{(e rr -e w f + 4e 2 r }] (5) 



and integrating / over the entire space, both within and 
outside the disk (to account for the deformation of the 
surrounding medium) and using the strains calculated 
from Eqs. 3 and 4, we get the energy E = 2ir(K + fi)Ar 2 
necessary to expand the disk by Ar. We may, now, elimi- 
nate the shape dependent prefactors by using the volume 
V = ttR 2 and the volume-change AV = 2-irRbAr of the 
disk, to obtain finally, the energy E — (K + fi)AV 2 /2V. 
Using the equipartition theorem we have therefore, 



{{AV b ) 2 )/V b = kT/[K(L b ) + n{L b )}, 



(6) 



relating the fluctuation of the volume V b = L\ of the sub- 
cell to the sum of the bulk and shear moduli. The above 
relation together with AV b /V b = e b xx (L b ) + e b yy (L b ) may 
now be used to obtain K{L b ) + fi(L b ) for our sub-cells. 

A similar treatment leads to a relation between \i and 
the local rotation of the system 9 (Eq.l). The rotation 
of a disk of radius R b by an angle 9 leads to an angular- 
displacement u v (r)=9-R 2 /r for r > Rb (u v (r) = 
for r < Rb). Applying Eq. 4 and integrating the energy 
density (Eq. 5) leads to the total energy for the rotation 
E = 2ix\i9 2 R\. Equipartition then yields, 
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V b ((29 b ) 2 ) 



(7) 



This equation has precisely the same structure as Eq. 6 
and, therefore, similar finite size scaling schemes can be 
applied. Thus Eq. 7 together with Eq. 6 enable the de- 
termination of the elastic constants of the system. 
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FIG. 3. (a)and (b)Probability distribution of the box-size 



averaged strain fluctuations of AVb/VJ, = e xx + e yy as a func- 
tion of the (linear) scaling box-size Lb = L/b for V = 508 (a) 
and 69 (b). (c) and (d) Same as in (a) and (b) for the angle 
@b = £ b xx ~ e yy for T = 508 (c) and 69 (d). Note the difference 
in the y-scale between (a) and (b) as well as between (c) and 
(d) 



In Fig. 3 the probability distribution of the incremental 
volume AVb/Vb — e b x + e b y and 9b are shown as a func- 
tion of the scaled box size b = L/Lb for two values of 
the interaction strength T. The width of the (Gaussian) 
distributions decreases both with increasing box-size Lb 
and increasing value of T. The mean square fluctuations 
are obtained by fitting the data to a normal distribution 
and extracting the standard deviation. 

Once the Lfc-dependent elastic moduli are obtained, 
usual finite size scaling [5] can be used to extrapolate 
the results to the thermodynamic limit. This is shown in 
Fig. 4 where we have plotted L b /L x l/[K(L b ) + n(L b )] 
as a function of L b /L. The slope of the curve, obtained 
from a straight line fit to the data, gives the value for 
1/[K(L — oo) + ji(L = oo)]. A similar finite size scaling 
is used to obtain \x(L — oo) separately. These results are 
shown in Fig. 2 both for K + (x - as obtained from Eq.6 
- and for [i from Eq.7 as a function of the interaction 
strength T. Note the accuracy of the determination of 
the elastic moduli for our system which is facilitated by 
the fact that the interaction potential T of our system 
is precisely calibrated. The straight lines through our 
data in Fig. 2 are the results from & T = calculation of 
the elastic constants of a 2D triangular solid composed of 
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particles interacting with an inverse cubic potential [16]. 
For all inverse power potentials the T = limit is exact to 
lowest order [17]. Considering the deformation of a per- 
fect static triangular solid of particles interacting with a 
r~ 3 potential one obtains the relation /i = K/\Q and the 
numerical result K — 3.461 • T, where the numerical co- 
efficient is evaluated by performing a rapidly convergent 
lattice sum. The agreement is excellent (considering the 
fact that no fitting parameter is available) over a wide 
range of interaction strengths T down to values of 70 - 
the melting transition occurs at T = 60 [7]. 




dislocation free regions of the sample — which is pos- 
sible only if the system is sufficiently far away from a 
melting transition [7]. Suggestions [5,12] to circumvent 
this problem are, however, computationally difficult to 
implement. The study of elastic properties of paramag- 
netic colloids in the presence of obstacles and inclusions, 
as well as dynamical elastic response is an interesting di- 
rection for further research. This is particularly suited 
for our technique since it provides a local (and therefore 
precise) probe [4] for elastic properties. 
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FIG. 4. Finite size scaling be- 

havior of [K(Lb) + fi(Lb)]~ ■ Lb/L (top) and ^(Lb)" 1 as a 
function of Lb/L dependent on the interaction strength P. A 
straight line fit to the curves gives the infinite system values 
of the (inverse) moduli. 

Finally, a few words on the possible uncertainties in- 
volved in our determination of the elastic moduli seem 
to be in order. Firstly, we have neglected all fluctuations 
of the magnetic moment, both in amplitude and angle. 
Since a super-paramagnetic colloid particle is of macro- 
scopic dimensions compared to typical magnetic length 
scales this assumption seems to be justified. Secondly, 
we have assumed that the particles fluctuate on a flat, 
two -dimensional air-water interface. An estimate [6] of 
the out-of-plane fluctuations is given by the ratio of the 
gravitational length l g = kT/mg (where m is the mass 
of the particles and g is the acceleration due to gravity) 
and the interparticle spacing; this is typically 1 in 10 3 . 
Lastly, a possible limitation of this scheme, at least in its 
present form, is the requirement that the displacement 
field u(r) be analytic in order to obtain strains by taking 
derivatives. Therefore one needs to restrict analysis to 
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